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| We discuss the physical properties of four quasar jets imaged with the Chan- 

dra X-ray Observatory in the course of a survey for X-ray emission from radio jets 
(Marshall et al. 2005). These objects have sufficient counts to study their spa- 

\q [ tially resolved properties, even in the 5 ks survey observations. We have acquired 

Australia Telescope Compact Array data with resolution matching Chandra. We 

^ | have searched for optical emission with Magellan, with sub-arcsecond resolution. 

The radio to X-ray spectral energy distribution for most of the individual re- 
gions indicates against synchrotron radiation from a single-component electron 
spectrum. We therefore explore the consequences of assuming that the X-ray 
emission is the result of inverse Compton scattering on the cosmic microwave 

k^J ' background. If particles and magnetic fields are near minimum energy density 
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in the jet rest frames, then the emitting regions must be relativistically beamed, 
even at distances of order 500 kpc from the quasar. We estimate the magnetic 
field strengths, relativistic Doppler factors, and kinetic energy flux as a function 
of distance from the quasar core for two or three distinct regions along each jet. 
We develop, for the first time, estimates in the uncertainties in these parameters, 
recognizing that they are dominated by our assumptions in applying the stan- 
dard synchrotron minimum energy conditions. The kinetic power is comparable 
with, or exceeds, the quasar radiative luminosity, implying that the jets are a 
significant factor in the energetics of the accretion process powering the central 
black hole. The measured radiative efficiencies of the jets are of order 10~ 4 . 

Subject headings: (galaxies:) quasars: general, galaxies: jets, X-rays: galaxies 

1. INTRODUCTION 

Following the remarkable discovery of an X-ray luminous, 100 kpc scale jet in PKS 0637- 
752 (Schwartz et al. 2000; Chartas et al. 2000), we have embarked on a survey to investigate 
the occurrence and properties of such systems. Initial goals were to assess the frequency of 
detectable X-ray fluxes from radio-bright jets, to locate good targets for detailed imaging 
and spectral followup studies, and where possible to test models of the X-ray emission by 
measuring the broad-band, spatially resolved, spectral energy distributions (SED) of jets 
from the radio through the optical to the X-ray band (Marshall et al. 2005; Schwartz et al. 
2003a,b). With an X-ray detection of 12 jets out of the first set of 20 observed, (Marshall 
et al. 2005), the survey has been successful in meeting these objectives. For four of these 
jets we have 40 to 130 total X-ray counts in our 5 ks observations. This suffices to construct 
broad-band spectral energy distributions, from which we can estimate magnetic fields, par- 
ticle densities, Doppler beaming factors and kinetic fluxes for independent, spatially distinct 
emitting regions, using the models of synchrotron radiation and inverse Compton (IC) scat- 
tering on the cosmic microwave background (CMB), which we adopt below (Schwartz et al. 
2003a,b). Deeper Chandra as well as HST observations have been approved for all these 
sources (Perlman et al. (2004), and in preparation). 

In parallel work, Sambruna et al. (2002, 2004) have undertaken a survey of 17 jetted 
radio quasars with Chandra and HST, with 10 exhibiting at least one knot in the Chandra 
images. We will compare our results with theirs in Section 4. 

Our survey is described in Marshall et al. (2005) (hereafter Paper I). We selected objects 
from two parent samples for which radio maps had been obtained with ~ 1" - 2" resolution 
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at 1-10 GHz: a VLA sample (Dec > 0°) of flat spectrum quasars with core flux densities 
S5GHZ > 1 Jy (Murphy, Browne, & Perley 1993) and an ATCA survey of flat spectrum Parkes 
quasars (Dec < -20°) (Lovell 1997), with core flux densities S2.7GHZ > 0.34 Jy. Selection was 
then made for objects with radio jets which extend beyond 2" from the core. The survey 
is comprised of sources for which we either anticipated extended structures with significant 
X-ray fluxes, based on scaling the extended 5 GHz flux using the X-ray to radio flux ratio of 
PKS 0637-752, (subsample "A"), or else selected by the morphological criteria of a one-sided 
linear jet (subsample "B"). PKS 0920-397 and PKS 1202-262 meet the "A" criterion, and 
PKS 0208-512, PKS 1030-357, and PKS 1202-262 meet the "B" criterion (Paper I). The 
definition of our parent samples by flat radio core spectra tends to select powerful sources 
with one-sided relativistic pc scale jets beamed toward our line of sight. The 5 or 2.7 GHz 
selection frequency emphasizes the jet rather than lobe emission. 

Section 2 presents the X-ray data and defines distinct spatial regions for further analysis. 
Section 3 gives the broad band spectral energy distribution, and Section 4 the physical prop- 
erties deduced by assuming that IC/CMB produces the X-ray emission. Section 5 discusses 
implications of the IC/CMB mechanism, and mentions alternate emission mechanisms. In 
Appendix A we estimate the systematic uncertainties in the magnetic fields and Doppler 
factors which we derive, and in Appendix B we present the basis for our calculation of the 
flux of kinetic energy carried by the jets. 



2. OBSERVATIONS OF THE JETS 

We observed the four quasars listed in Table 1 for about 5 ks, using ACIS-S. Data for 
the jets as a whole and X-ray results for the quasar cores appear in Marshall et al. (2005). 
We generally used the 1/4 sub- array mode to minimize pileup of the quasar core, but used 
the 1/8 subarray for PKS 0208-512 due to its greater flux. We requested roll angles such 
that the radio jet projected at least 30° away from the readout streak. Paper I (figure la, 
lg, lh, lk), shows the overlay of the 8.64 ATCA GHz radio contours on the X-ray images. 
Figure 1 shows the X-ray images of the four jets, and defines the regions used for the joint 
X-ray/radio spatial analysis. The jet regions were manually placed on the X-ray images, and 
therefore involve a certain amount of subjectivity. The regions are labeled R with numbers 
increasing away from the quasar. They are intended primarily to be regions larger than the 
instrumental resolution, and distinct from the quasar core, so we could derive independent 
model parameters characteristic of distinct volumes within each jet. 
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Fig. 1. — ACIS-S X-ray counts, 0.5 to 7 keV, binned into 0'.'492 pixels. The regions marked 
Rl, R2, R3 are used for the SED plots and the analysis of spatially resolved parameters. 
The circles centered on the quasars have a 1'.'23 radius, and enclose approximately 95% of the 
flux. The faintest color is one count per pixel. See text for definition of the marked regions. 
The straight lines give projected distances in the plane of the sky, at the redshift of each 
object. All images are in celestial coordinates with N up and E left, and cover 21"x 15". 
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The definitions of the spatial regions included a minimum of 5 X-ray photons. Back- 
ground counts are expected to be less than 0.1 in any region, and are ignored. Uncertainties 
are dominated by the X-ray statistics and model assumptions, as we shall discuss in Ap- 
pendix A, therefore we do not believe the exact region definitions affect any of the present 
conclusions. However, we do note with a question mark in Table 1 three regions for which 
there are less than 10 X-ray counts. The reality of detection and spatial location of these, 
and therefore their association with the radio jet, must be regarded as less certain. 

Table 1 gives the observed X-ray and radio data for each region. To derive the 1-keV 
X-ray flux density, fikev, and 2-10 keV rest frame luminosity, L x , we assume a power-law 
energy index of a=0.7, where f u oc v~ a . Our model that the radio and X-rays arise from the 
same simple power-law population of electrons implies that this is also the index of the radio 
emission. 1 For 3 of the quasars in Table 1 this is consistent with the spectral index from 4.8 
to 8.6 GHz, cufj, as tabulated in column 9, 2 within errors whose effect will be discussed in 
Appendix A. For PKS 0208-512 the indices deviate markedly from a=0.7; nevertheless we 
perform the formal calculation using the same index for comparison with the other objects. 
For this source, we have contamination of the Rl 4.8 GHz flux density by the core, and the 
R2 and R3 regions may be a hotspot and extended lobe, and so require different modeling. 
Our deeper Chandra and 20 GHz ATCA observations of this source will address these issues 
in the future. 

The X-ray properties reported here differ slightly from those in Paper I, because that 
paper considered the complete jet region, while in this work we omit some counts outside the 
regions marked in Figure 1. The angular sizes of the radio and X-ray regions are calculated 
by subtracting in quadrature the l'/2 FWHM of the 8.6 GHz images given in Paper I from the 
dimensions of each region. The regions are all resolved in the 8.6 GHz beam. We consider 
the regions as cylinders of angular length, 9 h and diameter 9 d . We use a flat, accelerating 
cosmology, with H =71 km s _1 Mpc -1 , Q m = 0.27, and fi A = 0.73. 



1 Note that the GHz emitting electrons have Lorentz factors ss 30 times larger than the electrons respon- 
sible for the IC X-rays 

2 Except for PKS 1202-262, the 4.8 GHz data are at poorer resolution than the 8.6 GHz images, and this 
produces additional errors in af $ 
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Table 1. Observations of Quasars, and X-Ray Jet Regions 
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Note. — Region designations followed with a question mark contain fewer than 10 X-ray counts, so 
their physical association with the radio emission is not certain. 

a From the NED database, operated by JPL for NASA 

b Rest frame 2 -10 keV luminosity, in units of 10 45 ergs s" 1 , assuming isotropic radiation 
°Spectral index from 4.8 GHz to 8.6 GHz, defined as oc v~ a 

d Angular length, 9i, and diameter, 9d, of the regions in Figure 1, after subtracting the radio beam 
l'/2 FWHM in quadrature 

c Quasar core properties 
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3. SPECTRAL DISTRIBUTIONS 



Figure 2 plots the spectral energy distributions for each of the regions shown in Figure 1. 
We have obtained Magellan optical observations of the four quasars (Gelbord et al. (2003, 
2004), Gelbord and Marshall, in preparation 2005). In general no statistically significant 
optical emission is detected. An exception is R2 of PKS 0920-397; however, there are other 
faint optical objects in the field and we cannot rule out a chance superposition at present. 
We therefore treat all optical data as upper limits. 
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Fig. 2. — Spectral energy distributions for the distinct regions of each jet. The downward 
arrows are our upper limit to the optical flux from any region. We use a circle, square, and 
triangle, respectively, to indicate regions progressively further from the quasar. 



For PKS 0920-397, PKS 1030-357 and PKS 1202-262 our upper limit g' magnitudes 
do not allow the X-rays to be a simple power-law extrapolation of the radio synchrotron 
emission. For PKS 0208-512 our limit of 25th magnitude prevents such an extrapolation 
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for R2. The 4.8 to 8.6 GHz radio index for PKS 0208-512 Rl and R3 also argues against 
extrapolation of the synchrotron emission to the X-ray region. 

We will present the interpretation of all four objects in the context of inverse Compton 
emission from the cosmic microwave background. This has been argued to be the most 
plausible mechanism for the X-ray jet in many powerful radio-loud quasars (Tavecchio et al. 
2000; Celotti et al. 2001; Marshall et al. 2001; Sambruna et al. 2001, 2002, 2004; Harris & 
Krawczynski 2002; Siemiginowska et al. 2002, 2003a,b). We discuss the consequences and 
limitations of this assumption, and mention alternate emission mechanisms in section 5. 

4. PHYSICAL PARAMETERS 

We assume that the X-ray emission from each region arises from inverse Compton scat- 
tering by the same power-law population of electrons, 71(7) = n 7~ m electrons cm -3 per 
unit 7, emitting the radio-synchrotron radiation from that region. The ratio of synchrotron 
to Compton power is just the ratio of energy density of the magnetic field to the energy 
density of the target photons, assuming the latter are isotropic in the jet rest frame (Felten 
& Morrison 1966). In applying that formalism to the powerful X-ray jets, one typically can- 
not find a credible source of target photons if one also assumes that the magnetic field and 
relativistic electrons are nearly in energy equipartition, and are not in relativistic motion 
(e.g., Schwartz et al. 2000). Tavecchio et al. (2000) and Celotti et al. (2001) resolved this 
dilemma by exploiting the enhanced apparent CMB density for electrons moving with bulk 
relativistic velocity, f3c ~ c, with respect to the isotropic CMB frame (Dermer & Schlickeiser 
1994). In the frame of a jet moving with bulk Lorentz factor T, the CMB energy density 
will exceed the magnetic field energy density at redshifts 

z > max [(0.556^ G /r - 1), 0], (1) 

(Schwartz 2002), where B m g is the magnetic field in micro-Gauss (l/zG=0.1nT), and T=(l- 
(3 2 Y l/2 . Under such conditions the IC/CMB may become the dominant energy loss mecha- 
nism for the relativistic electrons. The energy density due to the radiation field of a quasar 
emitting an isotropic, bolometric luminosity 10 44 L 44 ergs s -1 , falls below the CMB energy 
density at a distance 4.7 y / (L 44 /(l + z) 2 kpc from the quasar. For PKS 0208-512 this dis- 
tance is 34 kpc, so the quasar radiation field may play a role, especially in producing 7-ray 
emission. We will consider this in more detail in connection with our deeper Chandra and 
HST observations (Perlman et al., in preparation). For the other three objects, the quasar 
field drops below the CMB energy density 7 to 14 kpc from the quasar, which is much smaller 
than the deprojected distances we derive for the X-ray emitting regions. 
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From Felten & Morrison (1966), following their approximations that the flux density 
at any frequency is produced by a 5-function at a characteristic mean electron 7, and that 
magnetic fields, particles, and target photons are isotropic in the emitting region, the ratio 
of synchrotron to IC/CMB emission, extrapolated to some common frequency, will be 

OWSc) « ^ ^ (2) 

where p = T 2 (l + z) A po is the mean energy density of the CMB at redshift z in a frame 
moving with Lorentz factor T, po — 4.19 xlO -13 ergs cm -3 is the local CMB energy density, 
and the apparent temperature of the CMB in the jet frame is T = (1 + z)TT , where the 
local CMB temperature is T =2.728°K (Fixsen et al. 1996). There are two independent 
parameters among the direction to our line of sight, 9, the bulk Lorentz factor, T, and the 
effective Doppler factor <5=[F(1 — /3cos(#))] _1 . Since the asymptotic value of 5 is 2r for 
large T when the jet is beamed exactly in our direction, i.e., cos(0)=l, we make the common 
assumption r = S, in the absence of any information on 9. We assume the spectral index 
a=(m-l)/2=0.7, and use the known CMB parameters (Fixsen et al. 1996). We know that 
the CMB photons are highly anisotropic in the jet frame, but the photons we observe will 
be forward scattered by electrons moving near to the line of sight, so we approximate that 
they see the mean energy density. With these conditions, equation 2 gives an estimate of 
the magnetic field: 

B^g ~ 4.23 x 10- n 5(l + 2) 2 ' 18 (S synch /S IC ) a588 , (3) 

where Sic is the inferred X-ray flux density at 1 keV for the assumed spectral index a— 0.7, 
and S sync h is the measured radio flux density at 8.6 GHz. Appendix A discusses the sensitivity 
of our results to the spectral index. As noted above, the regions of PKS 0208-512 are not 
consistent with such a spectral index in the range 4.8 to 8.6 GHz. 

Without considering relativistic beaming, we can apply the minimum energy conditions 
to estimate the magnetic field (Moffet 1975): 

Bl » 328.8( (1 + fcl)jLsynch ) 2 / 7 Gauss (4) 
<pV 

This formula assumes a uniform magnetic field, and isotropic particle distribution. (Other 
approaches to the minimum energy formulation will result in a small dependence of the 
exponent on a; e.g., Worrall & Birkinshaw 2005). We assume unit filling factor, 0=1, and 
a ratio of proton to electron energy density, ki=l. We take the lower and upper limits of 
the observed radio spectrum to be z/i=10 6 Hz and z/ 2 =10 12 Hz, in order to integrate over v 
to get the total synchrotron luminosity, L synch . We consider the emitting volumes, V, to be 
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cylinders with the measured angular lengths 61 and diameters 64. Appendix A illustrates our 
sensitivity to these assumptions. In the jet frame we have B = B 1 /8 (Tavecchio et al. 2000; 
Harris & Krawczynski 2002; Dermer & Atoyan 2004), so that we have 

B5 » A00.8(L synch /V) 2 / 7 . (5) 

In Figure 3 we plot equations 3 and 5. The intersection of these curves for each feature 
in the four jets gives a solution for the unknowns B and 5. These values are tabulated in 
columns 3 and 4 of Table 2. We find magnetic fields of order 10 to 25 fiG and Doppler 
factors in the range 5 to 15 for the various regions in the jets. In Appendix A we quantify 
the uncertainties on those numbers. Compared to the analysis in Paper I, which considered 
the entire length of the jet and which assumed T = 10, we here use smaller volumes which 
naturally lead to larger values of the magnetic field Bi and smaller angles of the jet from our 
line of sight. 

Sambruna et al. (2004) have carried out a joint Chandra and HST survey, and report 
the analysis of ten knots which have been detected in both the optical and X-ray bands from 
six of their sources. They interpret nine of these in terms of the IC/CMB model, and derive 
magnetic fields in the range 3 to 12 fiG, and Doppler factors in the range 6 to 14. The values 
derived here are therefore quite similar, despite some differences in the model assumptions; 
e.g., the assumed shapes of the emitting region and values of r y m in- As noted by Tavecchio et 
al. (2004), the values of the Doppler factor are expected to be relatively robust since other 
derived quantities depend on powers of 5. 

From the magnetic field strength and our arbitrary assumption of a low-frequency cutoff 
at an observed v = 10 6 Hz, we can calculate the low energy electron spectrum cutoff 7 min 
which roughly corresponds to that frequency: ^ min = y/l(l+z)/(4.25B)). With that 7 miri we 
can calculate the total number density, n e , of relativistic electrons by equating the particle 
energy density with that of the magnetic field, which is approximately the minimum energy 
condition: 

(1 + ki) J^ max n m e c 2 7 1-m <i7 ~ B 2 /(8n). Columns 5 and 6 give r y m in and n e , respectively. 

For a fixed Doppler factor 5, the maximum angle by which the jet can deviate from 
our line of sight is arccosfy^ 2 — l)/5], which is also the angle for which 5=T as we have 
assumed. From this maximum angle, given in column 7 and the measured angular projection 
of the jet on the sky, we can compute the minimum intrinsic length of each jet region, as 
given in column 8 of Table 2. 
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Fig. 3. — The Tavecchio et al. (2000) style diagrams for the regions defined in Figure 1. The 
intersections of equations 3 (lines increasing with B) and 5 (lines decreasing with B) give 
solutions for the Doppler factor 5 and the magnetic field in the jet rest frame. The regions 
in order of increasing distance from the quasar are shown as dotted ("Rl"), dashed ("R2"), 
and solid ("R3"). 
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Table 2. Properties of the X-Ray Jets 
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a X-ray flux in Jet divided by X-ray flux in quasar 

b Calculated from B so that electrons of •y m in give 1 MHz synchrotron emission 
Calculated assuming bulk Lorentz factor T equals Doppler factor 5 
d Minimum distance from quasar, deprojected by l/sm(9 max ), in kpc 
c Kinetic power of jet, 10 46 erg s _1 
f De-beamed luminosity divided by kinetic flux, in 10~ 4 
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The kinetic energy flux carried by the jet in our observer frame is given by 

P jet = AT 2 (3c(w-p c 2 /r), (6) 

where A is the cross sectional area, T is assumed equal to 5, po is the rest mass density, 
and w is the total relativistic enthalpy density in the jet rest frame (e.g., Bicknell 1994). 
Calculation of this quantity is discussed in Appendix B, and it is tabulated in column 9 of 
Table 2. It gives the ability of the jet to do work on its surroundings, and explicitly excludes 
the rest mass energy of the particles, which does not normally enter the energy budget of 
the black hole and cannot normally be recovered from the jet. It differs by a correction of 
order unity from the commonly used bulk power calculated using only the internal energy 
density flux (e.g., Ghisellini & Celotti 2001). 



cm 
o 



o 



m 
O 



3 

X 





■ > 

CD 



10 :a _ 



A 



. PKS 0208 



1 =E 

_PKS 1030 



oo 



□ 



□ 



□ 



" PKS 1202 

<> 

- PKS 0920 



0.1 ' ' ' — ' 

1 00 1000 

Minimum Jet Distance 



[kpc] 



Fig. 4. — Kinetic power flowing through each region, plotted against the minimum depro- 
jected distance from the quasar. Larger symbols indicate regions of the jets in PKS 0208-512 
(triangles), PKS0920-397 (diamonds), PKS 1030-357 (squares), and PKS 1202-262 (circles). 
The dashed lines give the estimated radiative luminosity of the core quasar, with the symbol 
and name plotted to the left. 
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The kinetic energy flux carried by the jet in our observer frame is given by 

P jet = AT 2 (3c(w-p c 2 /r), (7) 

where A is the cross sectional area, T is assumed equal to 5, po is the rest mass density, 
and w is the total relativistic enthalpy density in the jet rest frame (e.g., Bicknell 1994). 
Calculation of this quantity is discussed in Appendix B, and it is tabulated in column 9 of 
Table 2. It gives the ability of the jet to do work on its surroundings, and explicitly excludes 
the rest mass energy of the particles, which does not normally enter the energy budget of 
the black hole and cannot normally be recovered from the jet. It differs by a correction of 
order unity from the commonly used bulk power calculated using only the internal energy 
density flux (e.g., Ghisellini & Celotti 2001). 

For a pure electron/positron jet, (i.e., ki—0), the powers would be a factor m 3 lower. 
For the case where charge neutrality was maintained by an equal number density of cold 
protons and with k\ ~ 10, the kinetic power would be about 20 times larger. We compare 
with the quasar bolometric radiative luminosity, estimated by fitting the radio loud template 
of Elvis et al. (1994) to the optical magnitude from the NED database, assuming isotropic 
emission, and integrating over all wavelengths. As shown in Figure 4, the kinetic powers are 
typically of order or larger than the bolometric accretion luminosity of the quasar. If the 
core optical emission is beamed, then the intrinsic luminosity of the cores are even smaller, 
and this conclusion is strengthened. As pointed out by Meier (2003), this requires that the 
jet formation be considered as an essential feature of the accretion process which is powering 
the quasar. The low efficiency, Table 2, column 10, with which these jets radiate their kinetic 
power is consistent with the ability to transport energy from the black hole core to distant 
radio lobes. 

For PKS 0208-512 Maraschi & Tavecchio (2003) have estimated a jet power pa 5xl0 47 
ergs s _1 at a few pc from the core. If we assume equal numbers of electrons and of cold 
protons, as they do, our average kinetic power for the regions Rl and R2 would be about 
1.2xl0 48 ergs s _1 . This is therefore consistent with the conclusion of Tavecchio et al. (2004) 
who find that the power transported to the outer jets of PKS 1510-089 and 1641+399 is 
similar to the energy flux traveling through the pc scale regions. 

5. DISCUSSION 

The kinetic flux carried by the jet equals the core luminosity for PKS 0208-512, and 
dominates for the other three quasars. From the discussion in Appendix B, we have for the 
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kinetic energy flux 

P jet » 1.25 x W 9 6 2 B 2 (6 d D A ) 2 (2 + x (6 - 1)/S), (8) 

where Da is the angular size distance and x is the ratio of rest mass energy density to 
particle enthalpy density, and we have considered a tangled magnetic field with B\ = \B 2 . 
The value of Pj et is not very sensitive to the values of B and 5 as long as f3 ~1, and we are 
near equipartition (with e c i cctr0 ns(l+ki)=B 2 /87r), since the product B5 is determined by the 
observed and assumed quantities in equation 4. The importance of the X-ray observations 
is that they show the jets must be in substantial relativistic bulk motion with (5 ~1, if we 
explain the X-rays as due to the IC/CMB mechanism. 

Alternate models for the X-ray emission have recently been considered. These are 
motivated at least in part by morphological arguments. Because the cooling length of the IC 
X-ray emitting electrons (7 ~ hundreds) is longer than that of the radio emitting electrons 
(7 ~ thousands), the extent of the X-ray emission should always be longer than that of the 
radio, in a model where the electrons are all accelerated at some discrete locations within a 
jet. However, some cases clearly show an X-ray peak upstream of a radio peak, which could 
be taken as evidence that the X-rays are produced by higher energy electrons; e.g., Stawarz 
et al. (2004). 

Dermer & Atoyan (2002) proposed that the observed X-rays in knots are due to syn- 
chrotron emission from electrons cooling by IC/CMB in the Klein-Nishina regime. This 
results in a high-energy (7 > 10 8 ) "hump" in the electron distribution function, manifested 
as a hardening of the synchrotron spectrum between UV and X-ray energies for an appropri- 
ately low magnetic field. In this model, the extrapolation of the harder X-ray spectrum to 
optical frequencies must always lie below the actual optical flux density. This model is ruled 
out by optical detections or upper limits at several knots; e.g., PKS 0637-752, (Chartas et 
al. 2000); Knot A of 1354+195, (Sambruna et al. 2004); Knot B of 1150-089, (Sambruna et 
al. 2002); and Knot C4 of 0827+243, (Jorstad & Marscher 2004). The limited statistics of 
our current X-ray data do not provide us with spectral information, and therefore we cannot 
check the validity of this model for the present quasars. 

It has also been suggested (e.g., Atoyan & Dermer 2004) that the X-rays are due to 
synchrotron radiation by a second, high-energy electron population. This requires a low- 
energy cut-off at a sufficiently high energy in the second electron distribution, so that its 
synchrotron emission does not over-produce the knot optical fluxes. Even if such a cut-off can 
be produced, the electrons will cool to energies below it in less than the escape time from the 
knot, to produce, in this regime, an electron distribution N e {^) oc 7~ 2 , resulting in a v~ x l 2 
optical spectrum that in many cases, such as in PKS 0637-752, overproduces the observed 
optical fluxes. In addition, given typical X-ray spectral indices (a x ~ 0.5—0.8; e.g., Sambruna 
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et al. 2004), and the fact that the cooling time of the X-ray emitting electrons is faster that 
the knot escape time, the injected electron distribution must have an index p flatter by one 
unit than the steady-state knot electron distribution index m, i.e., p = m — 1 = 2a x ~ 1 — 1.6, 
significantly flatter than that predicted by particle acceleration theories (p ~ 2 — 2.3, e.g. 
Kirk et al. 2000), unless there is much in situ acceleration in the knot. Spatially distributed, 
continuous acceleration (e.g., Jester et al. 2001, 2005; Marshall et al. 2002; Stawarz et al. 
2004; Perlman & Wilson 2005) may provide a solution to this problem . Such models can 
produce a pile-up of high energy electrons at the upper end of the electron distribution which 
could lead to X-ray synchrotron consistent with observations. However, similar to the model 
of Dermer & Atoyan (2002), such mechanisms cannot successfully reproduce emission by 
sources in which the extrapolation of the X-ray spectrum to optical frequencies lies above 
the observed optical flux. 

A different type of two-component model has been proposed by Aharonian (2002). 
According to this model, the X-rays are the synchrotron radiation of ultra high energy 
protons (7 P > 10 9 ). This requires a magnetic field B > 1 mG, larger by a factor of ~ 50 
than those typically used in leptonic models, and that the propagation of protons in the knot 
environment is taking place close to the Bohm diffusion limit. According to this model, the 
radio-optical component is synchrotron emission from an electron population injected in the 
knot, at a level which must be carefully fine-tuned to reproduce the radio-optical continuum. 

Although our objections to the above alternate explanations may not be decisive for 
these four quasars, the IC/CMB approach should be considered the simplest explanation as it 
adds only one parameter to the common assumption of minimum energy. The morphological 
data are certainly relevant and must be explained; however, they should not at present be 
considered as contradictory to the IC/CMB approach for two reasons. First, for the case that 
a radio peak appears downstream of an X-ray peak, we know that the naive model in which 
electrons are accelerated at a discrete location and then are advected downstream can not 
be directly applied, since in that case the peak emission must be coincident in both bands 
(Hardcastle et al. 2003), or at most be offset by an amount small compared to the angular 
resolution. Second, the cases where the X-ray emission peaks closer to the core, gradually 
decreasing outward, while the radio emission increases outward to peak practically at the 
jet terminus, might be explained if the large-scale jet gradually decelerates (Georganopoulos 
& Kazanas 2004) downstream from the first knot. The X-ray brightness then decreases 
along the jet because the CMB photon energy density in the flow frame decreases. At the 
same time, the deceleration leads to an increase of the magnetic field in the flow frame, 
which enhances the radio emission with distance. As a result the radio emission is shifted 
downstream of the X-rays and a rx increases along the jet. As shown in Figure lg of Paper 
I, PKS 0920-397 is similar to 3C 273 (Marshall et al. 2001) in displaying such morphology 
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and we note that our modeling (Table 2) of PKS 0920-397 gives a decreasing Doppler factor 
and increasing magnetic field going from Rl to R2, (but see Appendix A for systematic 
uncertainties). 

If the IC/CMB mechanism proves relevant to jet X-ray production, it may have the 
important cosmological implication that jets with identical intrinsic parameters will appear 
to have the same surface brightness at whatever redshift they might exist (Schwartz 2002). 
We note for the present jets, that even if they are not in relativistic motion and that their 
X-ray emission is not IC/CMB, in the equipartition models we have assumed the magnetic 
fields are in the range 50 to 100 /zG, (the product B5 from columns 3 and 4 of Table 2). For 
such magnetic fields, at redshifts larger than 3 to 4.5, respectively, the CMB will dominate 
the target photon energy density, and the predominant X-ray mechanism will be IC/CMB, 
providing the relativistic electron spectrum extends to sufficiently low Lorentz factors. The 
present jets would only be about 15 to 75 times fainter at such redshifts, so most would still 
be detectable in long but feasible Chandra observations of a few 100 ks. 

Improved radio imaging and spectra, as well as deeper X-ray observations, will allow a 
more quantitative confrontation of the alternate X-ray emission mechanisms. If we hypoth- 
esize that these jets are carrying a constant kinetic flux we expect the (B,S) to lie along a 
line of constant B5 if the average injection energy has been relatively constant. With radio 
spectral indices measured to ±0.05 we could distinguish 5's differing by about 2 (see Fig- 
ure 5). With similar constraints on the X-ray spectral indices, we could test the hypothesis 
that radio and X-ray arise from the same population of electrons. We hope to obtain such 
data in our upcoming ATCA and Chandra observations. 
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A. SYSTEMATIC UNCERTAINTIES IN DETERMINING B AND 5 

We consider the uncertainty in the values we derive for B and 5. We will apply equations 
3 and 5 to calculate alternate values of magnetic field for the case 5—1. We vary one 
parameter at a time to evaluate the range which each equation may give, and then transform 
the extremes by B oc 5 or B oc 1/5, respectively, to give a rectangular region representing 
the systematic uncertainty. This is intended to give a conservative error; however, since we 
cannot quantify the probability of the assumptions, we cannot assign a numerical confidence. 
We will illustrate by considering if we can distinguish the Doppler factors for the distinct 
regions Rl and R2 in PKS 0920-397. These are plotted as the solid squares in Figure 5. We 
will refer to this figure in the following discussion of systematic uncertainties. 

The largest uncertainly arises from the starting assumption of minimum energy. For 
lower magnetic fields the energy density U increases oc B~ 3 / 2 and for larger fields increases oc 
B 2 . For energy 10 times the minimum, the crosses give the resulting fields. If conditions are 
far from equipartition, the energy may be arbitrarily larger and we would have no constraints 
on the magnetic field. In such circumstances it is possible that the X-rays are produced by 
a different model; hence, we will not use this limit in constructing an uncertainly region for 
our assumed model conditions. 

In applying equation 4 we have used assumed values of the radio spectral index, the 
lower and upper radio spectrum cutoffs, the ratio of proton to electron energy densities, and 
the volume filling factor. The squares show the effect of varying the spectral index a from 0.6 
to 0.9. The open circles are for low frequency cutoffs of v x = 10 5 and 10 7 . Varying the upper 
frequency cutoff has an extremely small effect because the bulk of the energy is in the lowest 
energy electrons. The largest effect is due to varying the ratio of proton to electron energy 
density, which we assumed as equal, to the values k\= and 10, shown by the downward 
triangles. Similar errors would arise from a 50% error in the cylinder radius (which in turn is 
merely an assumed geometry). These define the extremes we will take for the uncertainty on 
the magnetic field estimated via the equipartition assumption. If the filling factor were an 
order of magnitude smaller, psO.I, B would be about twice as large, similar to the /ci = 10 
case. 

The uncertainty in the magnetic field estimated from the inverse Compton formalism is 
dominated by uncertainty in the radio spectral index. These are shown by the open squares 
at B ^0.5 and 3.3 for a = 0.5 and 0.9, respectively. The upward triangles show the effect of 
a factor of 2 uncertainty in the X-ray flux, due to the Poisson errors and the uncertainty of 
converting counts to X-ray flux density when we cannot measure the spectrum accurately. 
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Magnetic Field [fiG] 

Fig. 5. — An uncertainty analysis applied to the region Rl of PKS 0920-397. Solid dots 
give the magnetic fields according to equations 3 and 5 and for 5—1. Alternate values of 
the minimum energy magnetic field, B me , or the field to produce X-rays via IC/CMB, B IC , 
arise as follows: Crosses, internal energy 10 times larger than minimum; open circles, low 
frequency radio cutoff v\ = 10 5 or 10 7 ; downtriangles, ratio of proton energy to electron 
energy k\= or 10; open squares, electron energy index = 2.2 or 2.8; uptriangles, X-ray flux 
error of a factor of 2. The filled squares are the solutions given in Table 2 for regions Rl 
{B « 12) and R2 (B « 20). The enclosed rectangle about the point B = 12.1//G and 5=8.3 
gives our adopted region for the systematic uncertainty. 

For PKS 0920-397 we see that the net uncertainty in (B, 5) space about the solution for 
Rl almost extends to the R2 solution. Obviously the joint uncertainty regions would overlap. 
If we assume an independent structure for these two regions, then a single value of B and 
5 could be assigned. On the contrary, if the electron spectra are assumed to have the same 
shape and no other parameters varied greatly between the regions, then the uncertainties 
would be correlated and we would have evidence for deceleration of the jet. 
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B. CALCULATION OF THE KINETIC ENERGY FLUX 

There are various expressions in the literature for the Poynting flux and particle kinetic 
luminosities of relativistic jets. Some are incomplete and some are incorrect. One common 
misconception, for example, is that the energy flux density associated with the magnetic fields 
is equal to the Lorentz factor squared times the magnetic energy density times the velocity. 
Here we present the correct equations for jet power. Although these expressions result in 
factor of order unity corrections to commonly used expressions, it is useful to have precise 
expressions for these quantities to provide a common basis for comparisons of expressions 
for jet power. 

Following Bicknell (1994), let e be the total (rest-mass plus internal) energy density, 
Po the rest-mass density, w = e + p the relativistic enthalpy density, j3 — v/c the bulk jet 
velocity and T = (1 — /3 2 )~ 1//2 the bulk Lorentz factor. Then, in the most general case, the 
jet energy flux density, in the absence of magnetic fields, is given by 

F E>i = [wT 2 - r Po c 2 ] & (Bl) 

This first term is derived from the stress energy tensor of a relativistic ideal fluid. The second 
term subtracts the flux of rest-mass energy, which for the purposes of energy conversion in 
the external regions of AGN, is not generally relevant. 3 Note however, that the rest-mass 
energy contributes to the total energy flux, through its appearance in the the relativistic 
enthalpy, w. 

We consider the jet plasma to be a mixture of relativistically cold matter (p << e) and 
relativistic particles with pressure p = |e, where e = e — poc 2 is the internal energy density. 
The energy flux density becomes: 

(B2) 



r 2 cA 



r - i 



-Poc 



The synchrotron and inverse Compton radiation emitted by a jet depends upon the 
relativistic electron and positron energy density, so that we express the internal energy in 
terms of these quantities, and write: 

e = (l + fci)e e (B3) 
p c 2 = (l + k 2 )n e m e c 2 (B4) 



3 In a pure electron/positron jet the rest mass might be recovered via annihilation, and the last term in 
equation Bl should be dropped 
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where e e , n e and m e are the electron/positron energy density, number density and mass 
respectively and the parameters k\ and k 2 represent the contributions from ions to the 
internal energy density and rest-mass density. For an electron-positron jet in which only 
relativistic particles are present, k\ — k 2 — 0. For an electron-ion jet, hi, in principle could 
take a number of values, depending on the details of acceleration and re-acceleration and 



l + k 2 



fi 



1 + 



n e m e 



u 



(B5) 



where \i is the mean molecular weight, n» is the ionic number density for an ion of mass rrii 
and u is an atomic mass unit. For a plasma with normal cosmic abundances, 1 + k 2 ~ 2, 200. 



In terms of this parameterization, the energy flux density may be represented as 

1 



F E>i = T'cPi 
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-(l + fci)e e + — 



-(1 + k 2 )n e m e c 
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and, introducing the jet area A, the jet power (energy flux) attributed to the rest-energy 
density and internal energy density of the particles, is given by 



r 2 c pA 
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The rest-mass energy density of the relativistic electrons (n e m e c 2 ) may be determined 
from their energy density as follows. Suppose that the number density per unit Lorentz 
factor of the relativistic electrons are represented by a power-law in Lorentz factor, 7, for 

7min < 7 < 7 max ■ 

by 

"(7) = n ol~ a 



with a > 2. Then, 
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For 7 max > 7 min , n e m e c 2 /e e w (a - 2)(a - 1) 1 7 1 



min ' 



(B8) 

(B9) 
(BIO) 
(Bll) 



Note that the rest mass energy density of relativistic electrons and the internal electron 
energy density are only of the same order when 7 min ~ 1. However, for a jet of normal cosmic 
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composition, the rest-mass contribution to the jet power is comparable to or greater than 
the internal energy contribution for 7 min ;$ 1000. 

To calculate the Poynting flux we use a Cartesian coordinate system defined such that 
the x-axis is in the direction of the jet flow, the y-axis is perpendicular to the flow and the 
magnetic field is in the x — y plane. 

Let E and B be the electric and magnetic field in the lab-frame, with the corresponding 
fields E' and B' = B'^x. + B' ± y in the jet rest-frame. As a result of the high conductivity 
of the plasma, E' = 0. Therefore, in this coordinate system, the transformation from rest 
frame to lab frame (see Jackson (1975)) is given by 



E = T[3B' ± z 
and B = B',,Sc - FB' ± y 



(B12) 
(B13) 



The component of Poynting flux along the jet in cgs units is 
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For the case of a tangled magnetic field with uniform strength, then the average squared 
perpendicular magnetic field (B±) = 2/3{B 2 ) and the area — integrated jet electromagnetic 
power is given by: 

P5 M « AT 2 f3c ( ^) (B15) 



jet 



V 67T 



However, conservation of magnetic flux in a smooth expanding jet flow implies that B' ± will 
dominate over B[, in the outer jet. In that case we expect 
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In view of the above, the total jet power is given by: 

?/2 



P jet « T 2 c(3A 
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If the internal energy density and magnetic field are in equipartition, then (B' 2 )/8n 
(l + fci)e e . 
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